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Abstract. Heterogeneous, “aggregated” patterns in the 
spatial distributions of individuals are almost universal 
across living organisms, from bacteria to higher vertebrates. 
Whereas specific features of aggregations are often visually 
striking to human eyes, a heuristic analysis based on human 
vision is usually not sufficient to answer fundamental ques- 
tions about how and why organisms aggregate. What are the 
individual-level behavioral traits that give rise to these 
features? When qualitatively similar spatial patterns arise 
from purely physical mechanisms, are these patterns in 
organisms biologically significant, or are they simply epi- 
phenomena that are Jikely characteristics of uny set of 
Interacting autonomous individuals? If specitic features of 
spatial aggregations do confer advantages or disadvantages 
in the fitness of group members, how has evolution operated 
to shape individual behavior in balancing costs and benefits 
at the individual and group levels? Mathematical models of 
social behaviors such as schooling in fishes provide a prom- 
Ising avenue to address some of these questions. However, 
the Inerature on schooling models has lacked a common 
framework to objectively and quantitatively characterize 
relationships between individual-level behaviors and group- 
level patterns. In this paper, we brielly survey similarities 
and differences in behavioral algorithms and aggregation 
statistics among existing schooling models. We present 
preliminary results of our efforts to develop a modeling 
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framework that synthesizes much of this previous work, and 
to identify relationships between behavioral parameters and 
group-level statistics. 


Types of Aggregations 


Aggregation occurs in the smallest organisms—bacte- 
ria—and the largest—whales—and spans virtually the en- 
tire extant diversity of taxon, habitat, trophic level, life- 
history strategy, degree of mobility, and many other 
biological characteristics (Parrish and Edelstein-Keshet, 
1999; Camazine er a/., 2001). Physical aggregation can he 
regarded as part of a continuum in group integration. At one 
end of this continuum are territorial animals with little need 
to engage in information transfer and no need for group 
structure. At the other end are highly integrated, long-term 
associations between individuals that know—and_ perhaps 
are even related to—other members of the group, and in 
which members can potentially have high rates of direct and 
indirect information exchange. Honeybee hives, cetacean 
pods, and human communities are examples of these highly 
integrated groups (Wilson, 1975). In these systems, estab- 
lished pathways of long-term communication between 
known individuals (clones, siblings. reciprocating group 
members), or at predetermined locations (the hive, the calv- 
ing grounds, the dinner table), may supplement unmediate 
sensory contact. Thus, the members remain part of the 
“croup” even while they range widely in space. In these 
dispersed groups, coordinated function can be maintained as 
long as the necessary information is transferred between 
group members (e.g., the sensory integration systems of 
Norris and Schilt, 1988: Schilt and Norris, 1997), though 
distance between interacting components Inevitably affects 
the evolution and stability of emergent group properties 
(Hillier and Hanson, 1990; Latané er al.. 1995). 
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Between these extremes of group integration are what 
could be considered “prototypical” animal aggregations— 
herds, swarms, flocks, and schools. Within the fishes, over 
50% of species school—that is, display synchronous and 
coordinated movement—at some point in their life histories 
(Shaw, 1978), and an unknown additional number ageregate 
more coarsely. Prototypical aggregations exhibit coordi- 
nated motion, but group members are generally unrelated 
and never develop lasting relationships (in the game-theo- 
retic, “tit-for-tat” sense) with other members. Many of these 
groups are extremely large (e.g., a school of a billion 
herring). Individuals in such groups interact with a neigh- 
borhood of other members, but those may represent a van- 
ishingly small fraction of the group as a whole. This sug- 
gests that mechanisms which maximize information transfer 
among individuals could be evolutionarily beneficial. One 
example of this 1s a repeated arrangement within the group, 
reminiscent of crystal lattices, in which individuals assume 
preferred positions and orientations relative to their neigh- 
bors. Such arrungements could, for example, maximize sen- 
sory contact between members in such as way as to reflect 
ambient conditions, and the organisms’ predominant sen- 
sory systems, morphology, ete. (Parrish, 1992; Parrish and 
Edelstein-Keshet, 2000). 

We can characterize possible behavioral adaptations in 
members of these groups on at Jeast two levels: (1) short- 
term reactions to modify position with respect to immediate 
neighbors; and (2) behavioral responses that do not neces- 
sarily improve position relative to immediate neighbors but 
that contribute to group-level characteristics that ultimately 
benefit the individual by benefiting the group. These group- 
level adaptations are among the most fascinating—and the 
most difficult to assess—aspects of animal aggregations. 





Pattern versus Funetion 


Human perception tends to recognize attributes of the 
whole: an even density profile, polarity, distinct edges, or 
specific shape. While it is tempting to assume that conspic- 
uous features of biological aggregations are somehow ben- 
eficial, the existence of qualitatively similar patterns that 
arise from physical phenomena shows that this need not be 
the case. For example, some shapes found in three-dimen- 
sional schools (e.g., torus and funnel) are echoed in two- 
dimensional insect configurations (e.g., wheel), suggesting 
that such shapes may be adaptive for group members (Fig. 
1). However, these patterns could be evolutionarily neutral, 
or even pathological. Virtually identical shapes can be 
found in a wide range of inanimate aggregations, from water 
vapor to planets, in which these patterns arise from simple 
abiotic interactions between individual components, in the 
absence of evolutionary dynamics. Key steps in understand- 
ing biological aggregations in nature must be to distinguish 
biologically relevant features from nonadaptive epiphenom- 


ena, and to more explicitly and mechanistically describe the 
links between individual behaviors and group pattern. 

Dynamic patterns und movement are necessary charac- 
teristics of many biological aggregations, and are perhaps 
better criteria to distinguish adaptive responses from epi- 
phenomena. For example. fish schools display complex 
emergent properties such as coordinated motion and di- 
rected activity. Compression, hourglass, vacuole, fountain, 
and flash expansion (Fig. 2; Pitcher and Parrish, 1993) are 
all maneuvers that minimize predatory risk only if all mem- 
bers perform them correctly (e.g., Parrish, 1989). Parabolic 
formations in tuna schools may allow cooperative hunting 
advantages (Partridge et al, 1983). These emergent group 
properties appear to have readily apparent biological inter- 
pretations, and physical analogs may be harder to find. 
However, these dynamic group properties clearly confer 
more evolutionary advantages under some circumstances 
than others (e.g., evading small predators that target indi- 
viduals while becoming targets for large predators that 
target groups; Parrish, 1993). Furthermore, the key evolu- 
tionary question remains: do these behaviors involve trade- 
offs between short-term gain for the individual and long- 
term functioning of the group? If so, what is the 
evolutionary context in which selection for these behaviors 
occurs? 


Traffie Rules 


Assuming structure is advantageous, how is it main- 
tained? Laboratory and field attempts to address this ques- 
tion in fish schools have been limited (Partridge and Pitcher, 
1980; Aoki et al., 1986; Parrish and Turchin, 1997). in part 
because obtaining three-dimensional trajectories on specific 
individuals for a relevant period of time is difficult. Data 
that do exist are typically from highly artificial conditions 
(e.g., relatively small schools in highly lit still-water tanks; 
Parrish and Turchin, 1997). Three-dimensional tracking 
techniques have not yet advanced to the stage where it is 
feasible to observe large schools (i.e., over 10), In three 
dimensions, over long times (/.e., for more than seconds). 
Despite these difficulties, quantitatively accurate observa- 
tions of fish behavior within schools will undoubtedly be- 
come available in the next few years. However, while those 
data will provide a means of assessing short-term behavioral 
responses by fish to neighbors within a group, they will not 
by themselves provide the strong linkage between individ- 
ual and group characteristics that we require to understand 
the mechanics and evolution of schooling behaviors. 

Making this linkage requires an additional approach, 
namely, mathematical or computational models of school- 
ing behavior. These models posit a specific, quantitative set 
of behavioral interactions— essentially, they create a set of 
traffic rules—and quantitatively assess the emergent prop- 
erties of the resulting schools. Ideally, both the inputs (i.e., 
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Figure 1. Three- and two-dimensional expressions of emergent structure in animate and inanimate aggre- 
galions—whorl patterns. Left: tornado structure in tish (top) and water vapor (bottom). Middle and right: whorl 
and toroid in planets. fish. ants. and water vapor. clockwise from top left. Top left reprinted with permission by 
FPG; bottom left courtesy of National Severe Storm Laboratory; middle top and bottom courtesy of National 
Aeronautic Space Administration; top right taken by Norbert Wu, www.norbertwu.com © 1999; botiom right 
taken by T. Schneirla and reprinted with permission by W. H. Freeman & Company. 


individuals’ responses to neighbors) and model output 
(group-łevel characteristics) can be compared to data from 
real aggregations. 

A key purpose of modeling is to distinguish behavioral 
cause from organizational elfect by studying the conse- 
quences of various hypothetical social interaction rules. 
Most simulation models of animal aggregations im the lit- 
erature assign a set of forces that uct on the speed and 
direction of each individual and are modulated in response 
to other individuals or the local environment. Typical force 
components include locomotory (e.g., biomechanical forces 
such as drag), aggregative (e.g., long-range attraction, short- 
range repulsion), arrayal (e.g., velocity matching), and ran- 
dom (e.g., individual stochasticity; Griinbaum and Okubo, 
1994). The detailed biomechantes of locomotory forces are 
usually not considered in fish schooling simulations. In- 
stead, most simulations simply associate behavioral move- 
ment “decisions” with the movements that result. However, 
there are exceptions, notably a few modeling studies that 


specifically address hydrodynamic interactions between 
members of a school (e.g., Weihs, 1973). 

Published modeling studies have for the most part ad- 
dressed a relatively restricted range of behavioral algo- 
rithms. Most analyses focus on variation in one of three 
categories: behavioral matching, positional preference, and 
numerical preference. 

Behavioral matching (aka allelomimesis: Deneubourg 
and Goss, 1989) occurs when the individual agents try to 
match their behavior with other nearby agents. Most often, 
behavioral matching is modeled by having each agent ex- 
plicitly match its ortentation with that of its nearby neigh- 
bors (e.g., a zone of parallel orientation: Aoki, 1982; Huth 
and Wissel, 1992: Dill et al.. 1997; Table 1, Fig. 3). Less 
often, fish match their speed, either to that of their compan- 
ions, or to some arbitrarily decided value (e.g., Romey, 
1996; Vabo and Nøttestad, 1997). 

Positional preference refers to each fish having a pre- 
ferred position relative to one or more of its companions. 
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GN Hourglass 


Examples of coordinated movement and directed activity, both emergent properties of fish 


schools, which are also commonly cited defense tactics against predatory attack. Fig. 12.8. from Pitcher and 
Parrish 1993. Reprinted with permission from the author and Kluwer Academic Publishers. 


Usually, positional preference is formulated us a preferred 
distance to one or more nearest neighbor(s). Variations on 
positional preference include assigned distance zones (e.2., 
repulsion, parallel orientation, attraction, searching) within 
which neighbors are treated equally (Huth and Wissel, 
1992, Stocker, 1999) or continuous distance weighting 
(Warburton and Lazarus, 199]; Reuter and Breckling, 1994; 
Romey, 1996). In some models, other positional parameters 
influence responses, such as bearing angle to neighbors, or 
estimated collision time (Dill et al., 1997). The biological 
underpinning 1s obvious: that individual group members do 
not collide, that groups do not dissolve, and that stragglers 
join. Taken together, behavioral matching and positional 
preference describe what a fish should do, e.g., move to- 
wards or away from neighbors, align with neighbors, search 
for neighbors (Fig. 3; RPOA dependence, Table |), and how 
much consideration it should give any neighbor in the 
perception field (neighbor scaling rule. Table 1). 

Numerical preference refers to the number of neighbors 
to which a fish pays attention, which we generically refer to 
as the rule size. Variations include an a priori value (e.g., 4: 
Warburton and Lazarus, 1991) or a conditional value (e.g., 
choose up to 4 in the nearest zone: Aoki, 1982: choose up 
to 4, front prioritized: Huth and Wissel, 1992). Many simply 
have each fish average over all other fish (e.g., all visible: 
Reuter and Breckling, 1994: Vabo and Nøttestad. 1997; all 
fish: Romey. 1996). 

In general, modeling studies of schooling have been 
limited in several important respects, which future studies 
should aim to improve upon. Most models have used a 


relatively small population (e.g., 8 to 20 fish; Table 1), 
despite the fact that several studies point out that small 
numbers of fish may produce artificial results (e.g., Romey, 
1996). Most models have operated in two dimensions. Al- 
though this may be justified on the basis of computational 
resources, generalizing from two to three dimensions may 
not be straightforward. Many models include some stochas- 
tic or chaotic elements: however, the degree of replication in 
modeling studies is generally lower than ideal to character- 
ize the frequency distributions of possible outcomes. 

Most simulation studies also assume that all individuals 
are identical. Romey (1996) has shown that the inclusion of 
a single fish with different traffic rules will alter schooling 
(measured as group speed and turning rate). Furthermore, 
the inclusion of multiple agents with one of two rule sets 
produces sorting (Romey, 1996), which may translate into 
emergent properties (e.g., horizontal bands versus extended 
columns of ungulates, Gueron er al., 1996). The possibility 
that variability within the group may not only affect emer- 
gent properties such as coordinated movement and spatial 
pattern, but may also itself be an emergent property, clearly 
deserves further consideration in modeling studies. 

The parameter space explored by each model is typically 
only a small subset of possible variations, which ideally 
would include variations in (at least): initial position and 
velocity: the strength and type of stochastic components: 
spatial distribution of repulsion, parallel orientation, and 
attraction; and degree of variation between individuals in a 
group (Table |). Because most studies are not consistent in 
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Table 1 


Summary of fish school simulation parameters and output variables. R/PO/A-sequential zones of repulsion, parallel orientation (variably used), and 
attraction where “zone” implies equal force within a proscribed area (see Figure 3). Direction matching implies the use of a zone of parallel 
orientation. All other definitions as in the text. Numbers in parentheses refer to footnotes. 








Population R/PO/A Direction 
Reference Size Starting Orientation Velocity Dependence Matching? Neighbor Scaling Rule Size 
Aoki (1982, 1984) ie Random position & Random Zone (1,2) Yes Position, weighted by  Upto4 
orientation within front priority 
bounded urea 
Warburton & Lazarus 2-9 Regular square lattice Not specified Linear distance No Constant 1-8 
(1991) 
Huth & Wissel 8 Random position & Random Zone (1) res Constant; Single Up to 4, front 
(1990, 1992, 1994) orientation in fixed choice (4) prioritized 
area 
Reuter & Breckling 10, 20. 30, Random position, Not clear Linear distance yes Distance-weighted Al} visible 
(1994) 40, 50 orientation & average (I/D; 1/D°) 


speed within 
bounded area 








Romey (1996) 2-10 Random position & Constant Linear distance No Constant All 
orientation within 
bounded area (5) 
Vabo & Nøttestad 900 Random position & Constant Discrete distance No Single choice (7) All visible 
(1997) orientation within 
bounded & fixed 
areas 
Stöcker (1999) 12, 64 Random position & Constant Zone (1) nes Constant Variable (8) 
orientation within 
bounded area 
Random Component Aggregation Indices Aggregation Indices Aggregation Indices 
Reference (distribution) Individual Group Population 
Aoki (1982, 1984) Gamma (speed); Net Individual Displacement Compactness (1/SD of x, y 
Normal! (direction: 3) coordinates ) 
Warburton & Lazarus Circular normal later-Individual Distance Group Center Displacement: 
(1991) (direction) Group Shape (max 1]D°)/area 
Huth & Wissel Gamma (speed); Nearest Neighbor Distance Polarization, Expanse (rms 
(1990, 1992, 1994) Normal (direction) (neighbors 1—3) distance to center) 
Reuter & Breckling Normal (direction, Nearest Neighbor Distance Polarization, Expanse {mean 
(1994) speed) distance to center) 
Romey (1996) Uniform (direction; 6) Group Speed; Group Turning 
Rate 
Vabo & Nøttestad Discrete probabilities of Instability (maximum Number of Schools 
(1997) choosing less than individual distance traveled) 
optimal direction 
Stécker (1999) Uniform (direction; 9) Polarization; School Direction Number of Schools; 


School Size 





Note that random components selected from a normal distribution are centered on calculated means as per Figure 3. 

(1) Distance-specific zones of repulsion, parallel orientation, attraction, and searching. A blind spot (30-60 deg.) behind. 

(2) Beyond the zone of attraction, but within the visual angle (270 deg), attraction based on linear distance. 

(3) Outside the sensed area, direction was chosen from a uniform circular distribution (.e., random walk). 

(4) Single Choice = choose only one neighbor with a weighted probability N1, 1/2N2, 1/4N3, 1/8N4. 

(5) Individuals can sense entire area (~ 100 BL). 

(6) Every 15 or 30 time steps. 

(7) Highest (optimal) attraction value direction (of eight 45 deg increments) chosen. 

(8) All neighbors within the parallel orientation zone considered equally. Only closest neighbor within attraction zone considered. 
(9) Only applied to single fish with no neighbors (/.e., random walk searching). 
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Figure 3. Caricatures of the atiraction-repulston functions used by 
various authors to simulate schooling in fish (see also Table 1). Repulsion 
is always negative on the y-axis, whereas attraction is posilive. Zones of no 
repulsion or attraction (e.g., parallel orientation) are at zero. The x-axis is 
distance from fish of interest, displayed in arbitrary units. (a) Aoki (1982). 
(b) Warburton and Lazarus (1991). (c) Huth and Wissel (1990, 1992, 
1994). (d) Curves showing the relative strengths of repulsion (left), parallel 
orientation (center), and attraction (right) in Reuter and Breckling (1994). 
Note that total social force always sums to 1.0 (e) Romey (1996). (f) The 
attraction strength used by Vabo and Nøttestad (1997). (g) Stocker (1999), 
(h) Our simulation. 


these aspects (among others). direct comparisons of their 
results are not possible. 

Finally. although most authors quantified school charac- 
teristics in their simulations with one or more index of 
aggregation, there is little or no consensus on whieh indices 
are most biologically relevant. Consequently, different stud- 
ies tend to cite different statisties, again making eompari- 
sons difficult. For example, it is diffieult to determine 
whether Romey’s (1996) model provides the same amount 
of polarization as Huth and Wissel’s (1992), because the 
former did not report that particular emergent property. A 
related but more subtle problem is that no single paper has 
simultaneously presented indices at the individual, group, 
and population levels (Table 1). In part this is because the 
distinction between group and population is meaningless in 
some models (/.e., all individuals always interact and there- 
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fore the entire population is always part of the same group, 
e.g., Aoki, 1982: Reuter and Breckling. 1994). Even when 
this is not the case, most studies use initial conditions or 
behaviors that predispose individuals to start and remain 
within a single large group. However, there are again nota- 
ble exceptions, for example, the cellular automata models of 
Vabo and Nøttestad (1997) and Stocker (1999): see also 
Flierl ef al. (1999). 

Ineonsistencies between modeling approaches and aggre- 
gation indices have led different authors to draw seemingly 
contradictory conclusions about how individual behaviors 
affect group-level charaeteristies. For example, Aoki (1982; 
1984) reports that removing the repulsion zone causes 
school disintegration. whereas Huth and Wissel (1990, 
1992. 1994) assert it has little effeet on schooling. Huth and 
Wissel (1990, 1992, 1994) maintain that only a low number 
of influential neighbors (< = 3) are needed, whereas War- 
burton and Lazarus (1991) state that increasing the rule size 
is essential] for school stability. Reuter and Breckling (1994) 
and Romey (1996) found that adding random or individual 
movement had little effect on schooling, but Warburton and 
Lazarus (1991) concluded that the incorporation of random 
movements destroyed schooling. Because these studies lack 
a common framework for specifying behaviors and summa- 
rizing spatial patterns, it 1s not possible to resolve whether 
these reflect true hiological differences or simply compari- 
sons between apples and oranges. 


A Synthetic Simulation Study 


We are currently working on a simulation model of fish 
schooling that explores much of the parameter space cov- 
ered in previous models. Our model operates in three di- 
mensions, and does not (as yet) inelude any rules that, a 
priori, program emergent property (/.e., no behavior mateh- 
ing or zone of parallel orientation). Following Griinbaum 
and Okubo (1994) and others, we model the movement 
decisions of each fish as a sum of locomotory, social, and 
random forces. Model fish are based on our quantitative 
observations of schooling giant danios (Danio aequipinatus; 
unpub. data), with a mass of 1.7 g. a length of 5.3 cm, a 
maximum speed of 11.2 body lengths |BL]/s, and a maxi- 
mum aceeleration of 11.2 BL/s~). Locomotory foree eon- 
sists of drag, specitied as a drag coefficient (Cp) of 0, 0.01, 
or 0.02 (based on a range of estimated drag coefficients for 
streamlined bodies = 0.01-—0.001; Videler, 1993). Details 
of the methodology are presented elsewhere (Viseido et al., 
unpub. data). 

We have thus far used this model to study a range of 
population from 2 to 128 fish, but report here only results 
for 128 fish. We examined variation in three aspeets of the 
social foree: the attraction/repulsion funetion, the neighbor 
scaling rule. and the rule size. Attraction/repulsion was 
modeled as a linear, piecewise linear. or convex (sensu 
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Table 2 


Aggregation indices at the individual, group, or population level used to evaluate traffic rules in our sunulation 











Index Computation Level 
Curvature Mean turning angle per cm traveled Indiv dual 
Fractal Dimension How “dimensional” the path ts Individual 
NGDR Net displacement divided by grass displacement Individual 
NND Mean linear distance to nearest neighbor Individual 
Group Size Number ot fish within 5 body lengths of each other Group 
Expanse Mean quadratic distance to center for each group Group 
Polarization Mean angle devialion between individual and group heading Group 
Groups Number of groups in the universe, including stragglers Population 
Stragglers Number ot fish more than 5 body lengths from any other fish Population 
Collisions Number of collisions per fish per minute (reported as total collisions) Population 
Clark and Evans (1954) Index Observed NND vs. expected NND for a randomly spaced population Population 





Warburton and Lazarus, 1991) function with F = O cen- 
tered at 1.9 BL (Fig. 3h). We based this preterred nearest 
neighbor distance on schooling experiments with real fish 
from a previous study (Parrish and ‘Turchin, 1997). We used 
a perception distance of 180 BL (i.e... so individuals could 
detect neighbors to which they choose not to respond). A 
neighbor scaling rule weighted the attraction/repulsion force 
by the distance to each influential neighbor. Weighting was 
constant (Ze., equal weighting), negative linear, or negative 
sigmoid. Finally, the number of influential neighbors (rule 
size) was set at 4, 8, 16, or 24 fish. 

Random force was chosen for each timestep from a 
normal distribution with mean F = O and standard devia- 
tion equivalent to 1/3, 1/6, and 1/18 of the maximum social 
force an individual fish could experience from another fish. 
The amount of random force experienced by any individual 
was a function of the distribution of nearby neighbors. 
Individuals within a schoo] experienced relatively little so- 
cial force because most neighbors are at or near the equi- 
hbrium (ie. F = 0), thus random forces can be propor- 
tionately large. Individuals divorced from the group but 
close enough to sense it experienced maximum social force 
(i.e. attraction), and proportionately little random force. In 
our simulations, median random force as a function of total 
force ranged from 8%-18% depending on rule size. We 
computed fish movement as acceleration resulting trom the 
sum of random and social forces, minus drag. 

Each parameter combination in our simulations were 
replicated 15 times, for T800 time steps (ol ] s) each. We 
examined our model output with a suite of 1] aggregation 
indices that encompass individual- to population-level char- 
acteristics (Table 2). Statistics were based on time averages 
taken over the second half of each simulation (900 time 
steps), to avoid undue influence from initial conditions. We 
report here the results of systematic comparison within a 
single rule (e.g., drag off versus high drag) where the default 
conditions were medium random force, no drag, piecewise 
linear attraction/repulsion, constant neighbor sealing, and a 


rule size of 16. All numbers in text are means, and all 
statistics are one-way ANOVAs, 

In general, over the range of parameter values we used, 
asocial forces (drag, randomness) had a less drastic effect on 
fish schooling behavior than did the various functions com- 
posing the social force (Fig. +). Within the range of random 
force examined, there was no significant effect on any index 
of aggregation (Fig. 4). The addition of drag slowed fish 
down and increased the tortuousity of their paths, with a 
resulting increase in group size (without drag— 460 fish/ 
group, with drag— 104 tish/group) and cohesion (without 
drag—]}.5BL, with drag—t.1BL). These larger, slower 
groups milled rather than moving forward (net to gross 
displacement ratio [NGDR]| dropped trom 0.5 to 0.1). which 
decreased polarization. However. once added, there was no 
difference between the low and high drag condition (Fig. 4). 

Social forces had a significant effect on most measures of 
schooling. Fish with convex attraction-repulsion (A/R) re- 
sponses stayed closer together and had more tortuous paths 
than fish with piecewise linear A/R responses. The result 
was an increase in collisions (convex—720 collisions: 
piecewise linear—270 collisions; Fig. 4). Polarization was 
unalfeected: it was equally low (83 deg.) for both functions. 

Numerical preference had strong effects, both in terms of the 
number of influential neighbors and how that influence was 
scaled. Changing the neighbor scaling rule from equal weight- 
ing to a negative sigmotd Junction effectively increased the 
influence of close neighbors. Fish moved faster and in 
straighter paths. Schools became smaller (sigmoid—33 tish/ 
group; constant— 104 fish/group) and slightly more polarized 
(sigmoid—70 deg.; constant— 83 deg.), but also more diffuse 
(sigmoid—t.9BL; constant—1].1BL). As a result. collisions 
decreased (sigmoid—S collisions; coustant—270 collisions). 
Using a smaller rule size (4 fish) produced equivalent results to 
negative sigmoid neighbor scaling. Thus a higher number of 
influential neighbors increases school size at the expense of 
internal structure (measured as polarization). Paying attention 
to more neighbors leads to large groups because it causes a 
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Figure 4. Results of one-way ANOVAs systematically comparing 
parameter values within a given traffic rule listed as column headings. 
Within each traffic rule, arrows indicate a significant difference in the index 
between the listed condition (e.g., rule size, “4° nearest neighbors) to that 
of the default set in parentheses (one-way ANOVA, n = 15, P < 0.001). 


White diamonds indicate no significant difference. Rows are 11 mdices of 


aggregation. İndividual level: path curvature, fractal dimension (higher is 
more tortuous), net to gross displacement ratio (NGDR: higher ts 
straighter), nearest neighbor distance (NND). Group level: group size (in 
fish). mean distance to group center (or expanse, in body lengths), polar- 
ization (lower is more aligned). Population level: number of groups. 
number of collisions. number of stragglers, and the Clark and Evans (1954) 
index. See Table 2 for details on how each index was computed. 


“network of overlapping influences” (Warburton and Lazarus. 
1991): each group member “pulls on” a different set of com- 
panions. causing groups to remain cohesive. Over the range of 
rule size (RS) examined, 


Group Size = 3 RS — 7 (df = 59, =O. So) 


Several other higher order properties emerged during our 
modeling. There was a clear linear relationship between 
group speed and polarization, across all rule sizes. Faster 
groups are more polarized (Fig. 5a). We note that there ts no 
preferred speed, so that this 1s not an inevitable result of a 
polarized group, but probably reflects the fact that slower 
groups are composed of individuals with more tortuous 
paths. We also have begun to examine substructure within 
groups. At the level of an individual's immediate neighbors 
(here defined as N]J—6), nearest neighbor distance is virtu- 
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ally identical (Fig. 5b), indicating that at this scale, spacing 
is quite regular. However, when rule size is very large (24 
influential neighbors), this relationship breaks down (Fig. 
5b), perhaps indicating an upper limit for rule size. 
Although our results are preliminary, they highlight some 
interesting patterns. First. schoo! cohesion and coordinated 
movement can be partially achieved by positional and nu- 
merical preferences, without the need for behavioral match- 
ing. This shows that schools might indeed arise by large 
numbers of simple, overlapping interactions (Warburton 
and Lazarus, 1991; Romey, 1996). Our preliminary results 
also indicate that the shapes of the functions used to char- 
acterize both positional and numerical preferences are 
clearly very important (Fig. +). and worthy of careful in- 
vestigation in the future. Finally, our results demonstrate the 
need to consider many different quantitative indices of 
aggregation when assessing model results (Fig. 4). 


Emergent Evolution 


A central issue in modeling studies of aggregation is 
distinguishing between group- and population-level charac- 
teristics that are direct, inevitable consequences of the as- 
sumptions built into the model, and those that are truly 
“emergent properties.” In an evolutionary context, the im- 
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Figure 5. Higher-order index comparisons. (A) Group speed versus 
polarization for rule sizes 4, 8, 16, and 24. (B) Mean variance in the six 
nearest neighbor distances for individuals in the center of the school, over 
four rule sizes. Individuals were considered “central” if they were closer to 
the group center than the median distance lo center for the whole group. 
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portance of this distinction 1s that emergent properties may 
be closely tied to indtreet behavioral adaptations, which 
benefit individuals by benefiting the group. To satisfy the 
definition of emergent properties, the larger system (the 
group) must possess them, while its components (the indi- 
vidual members) do not (Clark et al., 1997). To a large 
extent, designating a group characterisuc as “emergent” or 
not is a matter of degree. For instance, models in which 
individuals actively align their directions to those of neigh- 
bors typically produce polarized groups (Huth and Wissel, 
1992 993: Reuter andi brechime, 1994. Stocker 11099). 
but this is probably too direct an outcome of the assump- 
tions to be considered an emergent property. Dill et al. 
(1997) pose a model in which each fish estimates its time to 
impact with other fish. In such a system, agents tplicitly 
consider their neighbors’ orientation, velocity, and position, 
rather than explicitly doing so. Collision avoidance would 
therefore seem to be a direct outcome and not an emergent 
property of this model. However, what about the converse: 
is polarization an emergent property of collision-avoidance 
behavior? 

From a biological perspective, un operational distinction 
might be to consider whether the group behavior benefits 
the members because of their membership. For example, 
Griunbaum (1998) used model results to suggest that when 
individuals simultaneously display gradient-climbing and 
alignment behaviors in a noisy environment, groups to 
which they belong more accurately climb environmental 
gradients, whereas loners can not. If displaying both behav- 
iors is costly to an individual, but it nonetheless benefits by 
being part of a group in a better environment, this would 
qualify as emergent under this criterion. 

Given the relatively underdeveloped state of schooling 
behavior models and observations that unify individual. 
group, and population characteristics, it is unavotdable that 
studies of the evolutionary “hows and whys” of schooling 
remain somewhat speculative and oversimplified. Nonethe- 
less, several studies have attempted to incorporate evolu- 
tionary ideas into schooling models by allowing iterative 
improvement of parameters, using an objective funetion as 
an analog to “fitness” in the evolutionary sense (Huse and 
Gjoesaete, 1999; Bonabeau et al., 2000). For instance, Hira- 
matsu ef al. (2000) link a standard stimulation of schooling 
(a la Huth and Wissel, 1992) to a genetic algorithm that 
selects the width of the parallel orientation zone, the size of 
the blind angle, the number of influential neighbors, and the 
neighbor sealing rule (front, side, or distance priority, ran- 
dom choice) to maximize fitness. In this case, fitness is 
modeled as the minimization of four output variables: near- 
est neighbor distance, polarization, expanse, and fractal 
dimension. As a result, parallel orientation zone ts maxi- 
mized and blind angle is minimized. A larger question is 
whether this definition of most fit—that is, least vari- 
able—is correct? 


Simulation studies like these hinge directly on two ques- 
tions (l) what ts optimal, and (2) what is the relationship 
between individual and group optimality? Such questions 
are certain not to have simple or universal answers. It is 
unlikely that we will soon attain complete, realistic models 
of an aggregating organism's life history that would allow 
us to truly represent “filness” in a population with evolving 
behaviors. However, we are optimistic about the role of 
schooling simulauons in addressing two central issues 
raised in this paper. First. we believe model studies will be 
successful in elucidating the mechanistic relationships be- 
tween individual-level behavior and the group-level spatial 
patterns they produce. Second, we expect that the combi- 
nation of explicit observations of schooling fish, improved 
simulations of the observed behaviors, and expanded use of 
statistical indices of aggregation on real and model aggre- 
gations will prove sufficient to identify the true behavioral 
algorithms used by fish. 
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